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Using a sample of 122 million T(3S') events recorded with the BABAR detector at the PEP-II 
asymmetric-energy e+e^ collider at SLAC, we search for the hi,{lP) spin-singlet partner of the 
P-wave Xbj{^P) states in the sequential decay T{3S) — > 7r°/if,(lP), hb{lP) 77?i,(15'). We observe 
an excess of events above background in the distribution of the recoil mass against the tt" at mass 
9902 ± 4(stat.) ±2(syst.) MeV/c^. The width of the observed signal is consistent with experimental 
resolution, and its significance is S.lcr, including systematic uncertainties. We obtain the value 
(4.3 ±1.1 (stat.) ±0.9 (syst.))xlO"'' for the product branching fraction B{r{3S) 7v"hb) x B{hb 
IVb)- 

PACS numbers: 13.20.Gd, 13.25.Gv, 14.40.Pq, 14.65. Fy 



To understand the spin dependence of qq potentials 
for heavy quarks, it is essential to measure the hyperfine 
mass splitting for P-wave states. In the non-relativistic 
approximation, the hyperfine splitting is proportional to 
the square of the wave function at the origin, vifhich is 
expected to be non-zero only for i = 0, where L is 
the orbital angular momentum of the qq system. For 
L = 1, the splitting between the spin-singlet ("'^-Pi) and 
the spin-averaged triplet state ((^Pj)) is expected to be 
AMhf = M((3Pj))-M(iPi) - 0. The ^Pi state of bot- 
tomonium, the hf,{lP), is the axial vector partner of the 



P-wave Xbj(l-P) states. Its expected mass, computed 
as the spin-weighted center of gravity of the Xbj{^P) 
states, is 9899.87± 0.27 MeV/c^ [ij. Higher-order cor- 
rections might cause a small deviation from this value, 
but a hyperfine splitting larger than 1 MeV/c^ might 
be indicative of a vector component in the confinement 
potential The hyperfine splitting for the charmo- 
nium ^Pi state he is measured by the BES and CLEO 
experiments d-lll to be -0.1 McV/c2. An even smaller 
splittingis expected for the much heavier bottomonium 
system [2|. 
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The hb{lP) state is expected to be produced in T{3S) 
decay via tt^ or di-pion emission, and to undergo a subse- 
quent El transition to the r]i,{lS), with branchingfrac- 
tion (BF) B{h{lP) IVbilS)) - (40 - 50)% [l, 0. 
The isospin-violating decay T{3S) TT^hb{lP) is ex- 
pected to have a BF of about 0.1% 0,0] j while theoret- 
ical predictions for the transition T(3S') — )- TT~^TT~hb{lP) 
range from ^ 10~^ up to ~ 10~^ 0. A search for 
the latter decay process in BaBAR data yielded an upper 
limit on the BF of 1.2 x lO"'' at 90% confidence level 
(C.L.) [13 . The CLEO experiment reported the 90% 
C.L. limit B(r(35) ^ T:°hb{lP)) < 0^27%, assuming the 
mass of the h to be 9900 McV/c^ 
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In this paper, we report evidence for the hi,{lP) state 
in the decay T{3S) — )■ 7r''/ib(lP). The data sample used 
was collected with the BaBAR detector [11] at the PEP- 
II asymmetric-energy e+e~ collider at SLAC, and corre- 
sponds to 28 fb~^ of integrated luminosity at a center- 
of-mass (CM) energy of 10.355 GeV, the mass of the 
T{3S) resonance. This sample contains (122 ± 1) mil- 
lion T(3i5) events. Detailed Monte Carlo (MC) simu- 
lations [13| of samples of exclusive T(35') Tr^hb{lP), 
hb(lP) 7776(15) decays (where the hb{lP) and 776(15") 
are hereafter referred to as the hb and the 77^), and of in- 
clusive T(35) decays, are used in this study. These sam- 
ples correspond to 34,000 signal and 215 million T{3S) 
events, respectively. In the inclusive T{3S) MC sample a 
BF of 0.1% is assumed for the decay r(35) ^ Tr° hb ^. 

The trajectories of charged particles are reconstructed 
using a combination of five layers of double-sided sili- 
con strip detectors and a 40-laycr drift chamber, both 
operating inside the 1.5-T magnetic field of a supercon- 
ducting solenoid. Photons are detected, and their ener- 
gies measured, with a CsI(Tl) electromagnetic calorime- 
ter (EMC), also located inside the solenoid. The BaBAR 
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detector is described in detail elsewhere 

The signal for T{3S) — ^ ir'^hb decays is extracted from 
a fit to the inclusive recoil mass distribution against the 
7r° candidates (77iiocoii(7r°)). It is expected to appear 
as a small excess centered near 9.9 GeV/c^ on top of 
the very large non-peaking background produced from 
continuum events (e+e^ — >■ qq with q = u,d,s,c) and 
bottomonium decays. The recoil mass, 77Zrecoii (''■'') = 
E*{ttO))^ -p*(7r0)2, where 



beam 



is the to- 



ViE, ______ 

tal beam CM energy, and E*{tt'^) and p*(7r") are the en- 
ergy and momentum of the 7r°, respectively, computed 
in the e+e~ CM frame (denoted by the asterisk). The 
search for an hb signal, requiring detection only of the re- 
coil tt", proved unfruitful because of the extremely large 
associated tt" background encountered. In order to re- 
duce this background significantly, we exploit the fact 
that the hb should decay about half of the time [3, |1| to 
77^6, and so require in addition the detection of a photon 
consistent with this decay. The precise measurement of 
the r]b mass [l3| defines a restricted energy range for a 
photon candidate compatible with this subsequent hb de- 



cay. The resulting decrease in hb signal efficiency is offset 
by reduction of the 7r° background by a factor of about 
twenty. A similar approach led to the observation by 
CLEO-c, and then by BES, of the he in the decay chain 
-0(25) —7- hcTT^ — >■ rjcJTT^ [sUg, where the 77c was identi- 
fied both exclusively (by reconstructing a large number 
of hadronic modes) and inclusively. 

The signal photon from hb — > 777b decay is monochro- 
matic in the hb rest-frame and is expected to peak at 
^490 MeV in the e+e" CM frame, with a small Doppler 
broadening that arises from the motion of the hb in that 
frame; the corresponding energy resolution is expected to 
be ^ 25 MeV. The Doppler broadening is negligible com- 
pared with the energy resolution. Figure 1 shows the re- 
constructed CM energy distribution of candidate photons 
in the region 250-1000 MeV for simulated T{3S) Tr°hb, 
hb -> 77^6 events before the application of selection cri- 
teria; the signal photon from hb jr]b decay appears as 
a peak on top of a smooth background. We select signal 
photon candidates with CM energy in the range 420-540 
MeV (indicated by the shaded region in Fig. 1). 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

E*(7) (GeV) 



FIG. 1: The reconstructed CM energy distribution of the can- 
didate photon in simulated T(3S) — > 7r°/ij,(lP), hb{lP) — > 
7771,(15) events. The shaded region indicates the selected 
E*{'y) signal region. 



We employ a simple set of selection criteria to sup- 
press backgrounds while retaining a high signal efficiency. 
These selection criteria are chosen by optimizing the ra- 
tio of the expected signal yield to the square root of the 
background. The T(3S') — > ir^hb, hb —5- 7776 MC signal 
sample is used in the optimization, while a small frac- 
tion (9%) of the total data sample is used to model the 
background. We estimate the background contribution 
in the signal region, defined by 9.85 < 77ii.ccoii ('i'^ ) < 9.95 
GeV/c^, using the sidebands of the expected hb signal 
region, 9.80 < 7nrocoii(7r°) < 9.85 GeV/c^ and 9.95 < 
TOrocoii(7r°) < 10.00 GcV/c2. 
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The decay of the rjb is expected to result in high final- 
state track multiplicity. Therefore, we select a hadronic 
event candidate by requiring that it have at least four 
charged-particle tracks and a ratio of the second to zeroth 
Fox- Wolfram moments 15 1 less than 0.6 [l6| . 

For a given event, we require that the well- 
reconstructed tracks yield a successful fit to a primary 
vertex within the e~^e~ collision region. We then con- 
strain the candidate photons in that event to originate 
from that vertex. 

A photon candidate is required to deposit a minimum 
energy in the laboratory frame of 50 MeV into a contigu- 
ous EMC crystal cluster that is isolated from all charged- 
particle tracks in that event. To ensure that the clus- 
ter shape is consistent with that for an electromagnetic 
shower, its lateral moment 17| is required to be less than 
0.6. 

A TT*^ candidate is reconstructed as a photon pair with 
invariant mass m{'jj) in the range 55-200 MeV/c^ (see 
Fig. 2). In the calculation of Wrccoii ('''''), the 7-pair in- 
variant mass is constrained to the nominal n'^ value [ij in 
order to improve the momentum resolution of the tt^. To 
suppress backgrounds due to misreconstructed 7r° candi- 
dates, we require |cos6'/i| < 0.7, where the hehcity angle 
9h is defined as the angle between the direction of a 7 



for each interval of ?Tirocoii(7r'^). 
xitf 



from 



candidate in the tt^ rest-frame, and the tt'^ 



direction in the laboratory. 

Photons from 7r° decays are a primary source of back- 
ground in the region of the signal photon line from 
hb Iflb transitions. A signal photon candidate is re- 
jected if, when combined with another photon in the 
event (72), the resulting 772 invariant mass is within 
15 MeV/c^ of the nominal 7r° mass; this is called a 7r° 
veto. Similarly, many misreconstructed 7r° candidates 
result from the pairing of photons from different tt^'s. 
A 7r° candidate is rejected if either of its daughter pho- 
tons satisfies the 7r° veto condition, with 72 not the other 
daughter photon. To maintain high signal efficiency, the 
TT*^ veto condition is imposed only if the energy of 72 in 
the laboratory frame is greater than 200 MeV (150 MeV) 
for the signal photon (for the 7r° daughters). With the 
application of these vetoes, and after all selection criteria 
have been imposed, the average tt*^ candidate multiplic- 
ity per event is 2.17 for the full range of 771(77), ^^id 1.34 
for the 7r° signal region (110 < 7^(77) < 150 MeV/c^). 
The average multiplicity for the signal photon is 1.02. 
For 98.4% of 7r° candidates there is only one associated 
photon candidate. 

We obtain the m^ccoi\{'^'^) distribution in 90 intervals of 
3 MeV/c^ from 9.73 to 10 GeV/c^. For each 771^0011 (tt^) 
interval, the 771(77) spectrum consists of a 7r° signal above 
combinatorial background (see Fig. 2). We construct the 
TOrecoii(7r°) spcctrum by extracting the 7r° signal yield 
in each interval of TMrccoiK"''^) from a fit to the 7*7,(77) 
distribution in that interval. The 777i.ecoii(7r'^) distribution 
is thus obtained as the fitted tt'^ yield and its uncertainty 
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FIG. 2: The result of the fit to the 771(77) distribution in 
data (data points) for the full range of mrocoii(7r''). The solid 
histogram shows the fit result, and is essentially indistinguish- 
able from the data; the shaded histogram corresponds to the 
background distribution. 



We use the MC background and MC 7r°-signal dis- 
tributions directly in fitting the 777(77) distributions in 
data For each 777recoii(7i'°) interval in MC, we obtain 
histograms in 0.1 MeV/c^ intervals of 777(77) correspond- 
ing to the 7r°-signal and background distributions. The 
Tr'^-signal distribution is obtained by requiring matching 
of the reconstructed to the generated tt^'s on a candidate- 
by-candidate basis (termed "truth-matching" in the fol- 
lowing discussion). The histogram representing back- 
ground is obtained by subtraction of the tt" si gnal from 
the total distribution. 

For both signal and background the qualitative 
changes in shape over the full range of rn-cccoi\{T^^) are 
quite well reproduced by the MC. However, the 7r° signal 
distribution in data is slightly broader than in MC, and 
is peaked at a slightly higher mass value. The 771(77) 
background shape also differs between data and MC. To 
address these differences, the MC 7r° signal is displaced 
in mass and smeared by a double Gaussian function with 
different mean and width values; the MC background dis- 
tribution is weighted according to a polynomial in 777(77). 
The signal-shape and background-weighting parameter 
values are obtained from a fit to the 777(77) distribu- 
tion in data for the full range of 7r7i.ocoii(7r°). At each 
step in the fitting procedure, the tt*^ signal and back- 
ground distributions are normalized to unit area, and 
a between a linear combination of these MC his- 
tograms and the 777(77) distribution in data is computed. 
The fit function provides an excellent description of the 
data (xV^-DF=1446/1433; 7Vi:)F=Number of Degrees 
of Freedom) and the fit result is essentially indistinguish- 
able from the data histogram. The background distribu- 
tion exhibits a small peak at the tt" mass, due to inter- 
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actions in the detector material of the type n7r+ — > pTr" 
or pTT~ — > riTT*^ that cannot be truth-matched. The nor- 
mahzation of this background to the non-peaking back- 
ground is obtained from the MC simulation, which in- 
corporates the results of detailed studies of interactions 
in the detector material performed using data [l^ . This 
peak is displaced and smeared as for the primary tt" sig- 
nal. 

The fits to the individual 771(77) distributions are per- 
formed with the smearing and weighting parameters fixed 
to the values obtained from the fit shown in Fig. 2. In 
this process, the MC signal and background distribu- 
tions for each mrecoiK''''^) interval are shifted, smeared, 
and weighted using the fixed parameter values, and then 
normalized to unit area. Thus, only the signal and back- 
ground yields are free parameters in each fit. The fit 
to the data then gives the value and the uncertainty of 
the number of tt" events in each mrocoii interval. The fits 
to the 90 771(77) distributions provide good descriptions 
of the data, with an average value of (x^/NDF) = 0.98 
(A''DF=1448), and r.m.s. deviation of 0.03 for the distri- 
bution of values. We verify that the fitted 7r° yield is con- 
sistent with the number of truth- matched tt^'s in MC to 
ensure that the 7r° selection efficiency is well-determined, 
and to check the validity of the 7r° signal-extraction pro- 
cedure. 

To search for an /if, signal, we perform a binned fit 
to the rrirecoiii'!^'^) distribution obtained in data. The hi, 
signal function is represented by the sum of two Crys- 
tal Ball functions [20| with parameter values, other than 
the hh mass, 771, and the normalization, determined from 
simulated signal T{3S) — > n'^hb events. The background 
is well represented with a fifth order polynomial function. 

Direct MC simulation fails to yield an adequate de- 
scription of the observed background distribution, al- 
though the overall shape is similar in data and MC. This 
is due primarily to the complete absence of experimen- 
tal information on the decay modes of the hb and rjb 
mesons. Simulation studies with a background compo- 
nent that is weighted to accurately model the distribu- 
tion in data show a negative bias of ~ 35% in the signal 
yield from a procedure in which the background shape 
and signal mass and yield are determined simultaneously 
in the fit. Consequently, we define a region of rrij-ecoai'^'^) 
chosen as the signal interval based on the expected mass 
value and signal resolution. The signal region includes 
any reasonable theoretical expectation for the hb mass. 
We fit the 777i.ocoii(7i'°) background distribution outside 
the signal interval and interpolate the background to the 
signal region to obtain an estimate of its uncertainty 
therein. Figure 3(a) shows the result of the fit to the 
distribution of m^cconi''^^) in data excluding the signal 
region, 9.87 < 777recoii(7r°) < 9.93 GeV/c^. The fit yields 
X^/NDF = 50.8/64, and the result is represented by the 
histogram in Fig. 3(a), including the interpolation to the 
hb signal region. 



We then perform a fit over the twenty intervals of the 
signal region to search for an hb signal of the expected 
shape. We take account of the correlated uncertainties re- 
lated to the polynomial interpolation procedure by creat- 
ing a 20x20 covariance matrix using the 6x6 covariance 
matrix which results from the polynomial fit. The error 
matrix for the signal region, E, is obtained by adding the 
diagonal 20 x 20 matrix of squared error values from the 
WrecoiKTr") distribution, and a value is defined by 

X^ - VE-W. (1) 

Here V is the column vector consisting of the difference 
between the measured value of the 77irGcoii(7r°) distribu- 
tion and the corresponding sum of the value of the back- 
ground polynomial and that of the hb signal function for 
each of the twenty 3 MeV/c^ intervals in the signal re- 
gion. In Fig. 3(b) we plot the difference between the 
distribution of TUrecoii (''■'') and the fitted histogram of 
Fig. 3(a) over the entire region from 9.73 GeV/c^ to 10.00 
GeV/c^; we have combined pairs of 3 MeV/c^ intervals 
from Fig. 3(a) for clarity. The yield obtained from the 
fit to the signal region is 10814±2813 events and the hb 
mass value obtained is m = 9902 ± 4 MeV/c^ with a x^ 
value of 14.7 for 18 degrees of freedom. 

In order to determine the statistical significance of the 
signal we repeat the fit with the hb mass fixed to the spin- 
weighted center of gravity of the Xbj{^P) states, 777 = 
9900 MeV/c^. The signal yield obtained from the fit is 
10721 ± 2806. The statistical significance of the signal, 
calculated from the square-root of the difference in x^ 
for this fit with and without a signal component is 3.8 
standard deviations, in good agreement with the signal 
size obtained. 

Fit validation studies were performed. No evidence 
of bias is observed in large MC samples with simulated 
hb mass at 9880, 9900, and 9920 McV/c^. In addition, 
the result of a scan performed in data as a function of 
the assumed hb mass indicates that the preferred peak 
position for the signal is at 9900 MeV/c^, in excellent 
agreement with the result of Fig. 3(b). 

We obtain an estimate of systematic uncertainty on the 
number of tt'^'s in each TrarocoiK''"*^) interval by repeating 
the fits to the individual 771(77) spectra with the line- 
shape parameters corresponding to Fig. 2 varied within 
their uncertainties. The distribution of the net uncer- 
tainty varies as a third order polynomial in 777i.ecoii(7''°). 
We estimate a systematic uncertainty of ±210 events on 
the hb signal yield due to the 7r°-yield extraction pro- 
cedure by evaluating this function at the fitted hb mass 
value. 

The dominant sources of systematic uncertainty on the 
measured hb yield are the order of the polynomial de- 
scribing the 7r7rocoii(7r°) background distribution, and the 
width of the hb signal region. By varying the polyno- 
mial from fifth- to seventh-order, and by expanding the 
region excluded from the fit in Fig. 3(a) from (9.87-9.93) 
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FIG. 3: (a) The mi.ocoii(7r°) distribution in the region 9.73 < 
77irccoii(7r*^) < 10 GeV / c? for data (points); the solid histogram 
represents the fit function described in the text. The data in 
the hb signal region have been excluded from the fit and the 
plot, (b) The mi.ccoii('?r'') spectrum after subtracting back- 
ground; in the hb signal region the data points are shown as 
squares, and the area with diagonal shading represents the 
uncertainties from the background fit; the shaded histogram 
represents the signal function resulting from the fit to the 
data. 



After combining these systematic uncertainty esti- 
mates in quadrature, we obtain an effective signal sig- 
nificance of 3.3 standard deviations. The smallest value 
of the significance among those calculated for the varied 
fits in the systematics study is 3.1 standard deviations. 
The yield is 10814 ± 2813 ± 1652 events and the hb 
mass value m = 9902 ± 4 ± 2 MeV/c^, where the first 
uncertainty is statistical and the second systematic. The 
resulting hyperfine splitting with respect to the center of 
gravity of the Xbj{^P) states is thus AA/hf = -I-2±4±2 
MeV/c^, which agrees within error with model predic- 
tions 0,11. 

To convert the /if, signal yield into a measurement of 
the product BF for the sequential decay T(35') — >■ 7r°/i;,, 
hb 777b, wc determine the efficiency £5 from MC by 
requiring that the signal tt" and the 7 be truth-matched. 
The resulting efficiency is ts = 15.8 ± 0.2%. Monte 
Carlo studies indicate that photons that arc not from an 
hb — > 7^76 transition can satisfy the selection criteria when 
only the T(3S') — > ir^hb transition is truth-matched. This 
causes a fictitious increase in the hb signal efficiency to 
e = 17. 9± 0.2%. Therefore, the efficiency for observed 
hb signal events that do not correspond to hb — ^ 77/6 
decay is Ae = 2.1%. However, there is no current experi- 
mental information on the production of such non-signal 
photons in hb and rjb decays. Furthermore, the above es- 
timate of efficiencies in MC docs not account for photons 
from hadronic hb decays, since the signal MC requires 
hb — > l^b- We thus assume that random photons from 
hadronic hb decays have the same probability Ae to sat- 
isfy the signal photon selection criteria as those from r/f, 
decays. We assume a 100% uncertainty on the value of 
Ae when estimating the systematic error on the product 
BF. 

Wc estimate the product BF for T(iS) ir^hb, 
hb — > ji]b by dividing the fitted signal yield, N, corrected 
for the estimated total reconstruction efficiency, by the 
number of T{3S) events, Nx{3S), in the data sample. We 
obtain the following expression for the product BF: 



B{T{3S) ^ n'^hb) X B{hb ^ 7^/6) 



N 



N- 



r(3S) 



1 



es C" 



(2) 



GeV/c^ to (9.85-9.95) GeV/c^, we obtain systematic un- 
certainties of ±1065 events and ±1263 events, respec- 
tively, taken from the full excursions of the hb yield un- 
der these changes. Similarly, wc obtain a total systematic 
uncertainty of ±1.5 McV/c^ on the hb mass due to the 
choice of background shape. 

The systematic uncertainty associated with the choice 
of signal lineshape is estimated by varying the signal 
function parameters, which were fixed in the fit, by ±1(T. 
We assign the largest deviation from the nominal fit re- 
sult as a systematic error. Systematic uncertainties of 
±154 events and ±0.3 McV/c^ arc obtained for the hb 
yield and mass, respectively. 



where 



C = 1 



Ae 



es B{hb jT]b) 



(3) 



is the factor that corrects the efficiency es for the non- 
signal hadronic hb and rjb contributions. In this equation, 
we assume a BF value B{hb —J- jrjb) = 45 ± 5% accord- 
ing to the current range of theoretical predictions. The 
corresponding correction factor is 1 — C '-^ 30%, with a 
systematic uncertainty dominated by the uncertainty on 
Ae. 

We obtain B{r{3S) n^hb) x B{hb Tib) = 
(4.3 ± 1.1 ± 0.9) X 10"^, where the first uncertainty is 
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statistical and the second systematic. The result is con- 
sistent with the prediction of Ref. which estimates 
4 X 10^^ for the product BF. Since the /ife-decay uncer- 
tainty reduces the significance of the product BF rela- 
tive to that of the /if, production, we may also quote 
an upper limit on the product BF. From an ensem- 
ble of simulated events using the measured product BF 
value, and the statistical and associated systematic un- 
certainties (assumed to be Gaussian) as input, we obtain 
B{T{SS) TT^hb) X B{hb jT]b) < 6.1 X 10-4 at 90% 
C.L. 

In summary, wc have found evidence for the decay 
T{3S) 7r°/ib, with a significance of at least 3.1 standard 
deviations, including systematic uncertainties. The mea- 
sured mass value, m = 9902±4(stat.)±2(syst.) MeV/c^, 
is consistent with the expectation for the hb{lP) bot- 
tomonium state 0, [2l| , the axial vector partner of the 
Xb.iilP) triplet of states. We obtain B{r{3S) ■n^hb) x 
B{hb 777fc) = (4.3 ± 1.1 (stat.)±0.9 (syst.))xlO-4 
(< 6.1 X 10-4 at 90% C.L.). 
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Note added in proof: After this paper was sub- 
mitted, preliminary results of a search for the hb in the 



reaction e^t 



hb{nP)TT^TT in data collected near 



the T{5S) resonance have been announced by the Belle 
GoUaboration (g^]. The hb{lP) mass measured therein 
agrees very well with the value reported in this paper. 



